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Abstract 

This work presents a first study of a radiative calculation for the gravita- 
tional axial anomaly in the massless Abelian Higgs model. The two loop 
contribution to the anomalous correlation function of one axial current and 
two energy- momentum tensors, (A a (z)T IIU (y)Tpt T (x)), is computed at an or- 
der that involves only internal matter fields. Conformal properties of massless 
field theories are used in order to perform the Feynman diagram calculations 
in the coordinate space representation. The two loop contribution is found 
not to vanish, due to the presence of two independent tensor structures in the 
anomalous correlator. 

PACS: 04.62. +v, ll.10.-z, 11.15.Bt, 11.25.Db, 11.25.Hf, 11.30.Rd, 11.40.Ha 
Keywords: Anomalies, Perturbation Theory, Conformal Symmetry, Quantum Gravity 



Typeset using REVTpjX 



*E-mail address: pachos@ctp.mit.edu 
^E-mail address: ricardos@ctp.mit.edu 



I. INTRODUCTION AND DISCUSSION 



Perturbation theory anomalies have been known for a long time, starting with the work 
by Bell and Jackiw [jl] and by Adler concerning anomalies in gauge theories. The fact 
that there are no radiative corrections to the one loop result for the anomaly has been 
countlessly proven or brought into question since the early work of Adler and Bardeen 
For this reason, explicit calculations of possible radiative corrections to the one loop anomaly 
are of particular interest. 

Using a method proposed by Baker and Johnson ||, Erlich and Freedman recently per- 
formed such an explicit calculation for the two loop contribution of the anomalous correlation 
function (A /I (x)A u (y)A p (z)) of three chiral currents, in the Abelian Higgs model and in the 
Standard Model ||. In here, we wish to extend such calculation to the case of a gravitational 
background. 

The Adler-Bell- Jackiw anomaly concerns the divergence of the axial current in a gauge 
field background. The calculation of the divergence of the axial current in a gravitational 
field background was later performed by Delbourgo and Salam |J, Eguchi and Freund |7j 
and Delbourgo ||. As in the ABJ case, these authors found an anomaly associated to the 
conservation of the axial current, the gravitational axial anomaly. Later, Alvarez-Gaume 
and Witten showed the significance of gravitational anomalies for a wide variety of physical 
applications ||. 

The question of absence of radiative corrections to the one loop result obtained for the 
gravitational axial anomaly is an issue not as well established as it is in the gauge theory 
case. This is the reason why we proceed to perform an explicit two loop calculation, adopting 
the spirit in |J. However, calculating the two loop contribution to the gravitational axial 
anomaly is a much longer task than to do so for the gauge axial anomaly. In this paper we 
shall address the first part of the computation, by calculating the abnormal parity part of 
the three point function involving one axial vector and two energy-momentum tensors at a 
specific two loop order in the Abelian Higgs model. The reason we choose to work in this 
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model is due to the recent interest arising from the gauge anomaly case in , and also due 
to the fact that this model is a simplified version of the Standard Model. In order to set 
notation, the anomalous correlator we shall be dealing with is: 

(AaWTMTpvW), (1-1) 

where A a is the axial current and T^ v the energy-momentum tensor. 

The method of calculation [|||| is based on conformal properties of massless field theories, 
and also involves ideas from the coordinate space method of differential regularization due 



to Freedman, Johnson and Latorre HlQ[ . In particular, the correlator (1.1) will be directly 
calculated in Euclidean position space and a change of variables suggested by the conformal 
properties of the correlator will be used in order to simplify the internal integrations. The 
order in two loops we shall be working involves no internal photons, but only internal 
matter fields (the scalar and spinor fields in the Abelian Higgs model). However, in this 
case diagrams containing vertex and self-energy corrections will require a regularization 
scale. To handle this technicality we shall introduce photons in our calculation, as there 
is a unique choice of gauge fixing parameter (in the photon propagator) which makes both 
the self-energies and vertex corrections finite. These "finite gauge photons" are merely a 
technical tool employed in the calculation. 

The use of conformal symmetry to construct three point functions is well established. Of 



particular interest to us is the work by Schreier jy], where three point functions invariant 
under conformal transformations were constructed. For the case of one axial and two vector 
currents, it was shown that there is a unique conformal tensor present in the three point 



function. More recently, Osborn and Petkos [|12|] and Erdmenger and Osborn [|T^] have 



used conformal invariance to compute several three point functions involving the energy- 
momentum tensor. However, the case of one axial current and two energy- momentum tensors 
was not considered. 

What we find in here is that, even though at one loop there is only one conformal tensor 
present in the correlator (1.1) - the one that leads to the contraction of the Riemann tensor 



with its dual in the expression for the anomaly -, at two loops there are two independent 
conformal tensors present in the correlator. This is unlike the gauge axial anomaly case 
where the only possible tensor is the one that leads to the field strength contracted with its 
dual in the anomaly equation. Precisely because of the presence of these two tensors in the 
two loop result for the three point function, this correlator does not vanish. Again, this is 
unlike the gauge axial anomaly case ||. 

The two linearly independent conformal tensors present in the anomalous correlator 
are the ones in expressions (3.6) and (3.7) below (where the notation is explained in the 
paragraphs leading up to these formulas). One thing we would like to stress is that every 
diagram relevant for our calculation is either a multiple of one of these tensors, or a linear 
combination of them both. 

Two comments are in order. First, the existence of two independent tensors in the two 
loop correlator could seem to indicate the existence of a radiative correction to the anomaly. 
On the other hand, the fact that the correlator does not vanish at two loops does not mean 
that its divergence (the anomaly) does not vanish at two loops. 

Another point of interest is to follow [ID| and study the differential regularization of 



the one loop triangle diagram associated to the gravitational axial anomaly, Figure 1(a). 
This is done in the Appendix. What one finds is that differential regulation entails the 
introduction of several different mass scales. Renormalization or symmetry conditions may 
then be used to determine the ratios of these mass scales. In the gauge axial anomaly case 



it was found that there is only one mass ratio [10]. In this gravitational axial anomaly case 



we have shown in the Appendix that there is more than one mass ratio. This multiplicity of 
the mass ratios introduces new parameters that could be able to cancel all potential (new) 
anomalies. Apart from presenting part of these different scales we shall not proceed with 
their study. Here, we shall only restrict to the calculation of the correlation function, which 
by itself consists a lengthy project. Extracting the two loop contribution to the gravitational 
axial anomaly from our three point function is a question we hope to report about in the 
near future. 



The structure of this paper is as follows. In section 2 we present the massless Abelian 
Higgs model, as well as a review of the basic ideas involved in the method of calculation we 
use. This includes the calculation of the one loop triangle diagram. Then, in section 3 we 
perform our two loop calculation, with emphasis on rigorous details. The many contributing 
diagrams are organized into separate groups, and then analyzed one at a time. 

II. THE ABELIAN HIGGS MODEL AND CONFORMAL SYMMETRY 

We shall start by presenting the massless Abelian Higgs model. In four dimensional 
Euclidean space, its action is given by: 

S = J d 4 x{^F, u F, u + (D^Drf + 4>i»D^ - f$(I4 + Rtf)il> - ^W}, (2.1) 

where we have used L = |(1 — 75) and R = |(1 + 75), with 75 = 71727374- The covariant 
derivatives are: 

D^cj) = (dp + igAJfa 

= (dp + \igA^ b )i), (2.2) 

so that the theory is parity conserving with pure axial gauge coupling. 

Next we introduce a background (external) gravitational field, in order to properly define 
the energy-momentum tensors associated to the scalar and spinor matter degrees of freedom. 
A simple way to do this is to couple our model to gravity, so that a spacetime metric g^ u {x) 
is naturally introduced in the Lagrangian as a field variable. Then we can obtain the 
energy-momentum tensor by varying the Lagrangian with respect to the metric g^x) as 
T^x) = ^j^Jx) I d 4 x y/^g C, where T^ ll/ (x) is manifestly symmetric. In addition we have 
to ensure that it is conserved and traceless, obtaining finally for the fermion field, 

72,=^70A)^, (2-3) 

and for the boson field, 

T* = \{d^ + <%0 f - U, v 8 a ^ d a <p - l -{4> d&tf + ^ d&t)}, (2.4) 
where (fxu) = \iv + v\i. 
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One should observe that in the two loop calculation we are interested in computing the 
order 0(gf 2 k 2 ) correction to the correlator, where g is the gauge coupling, / the scalar- 
spinor coupling, and k the gravitational coupling. This means that there are no internal 
photons in the associated diagrams, as these would be of order 0(g 3 k 2 ) - we shall only 
need photons as the external axial current, and in order to handle some of the potential 
divergences in the calculation (see section 3). This is why in (2.3) and (2.4) the scalar and 
spinor matter degrees of freedom are decoupled from the gauge field. 

Conformal symmetry plays a central role in our calculations, as it motivates a change 
of variables that simplifies the two loop integrations. Due to the absence of any scale, our 
model is conformal invariant. The conformal group of Euclidean field theory is 0(5, 1) ||12|| . 
All transformations which are continuously connected to the identity are obtained via a 
combination of rotations and translations with the basic conformal inversion, 



X>M x' 2 ' 



|f=x 2 (^-^)^x 2 ^(x). (2.5) 

The Jacobian tensor, J^ v (x), which is an improper orthogonal matrix satisfying J^x) = 
J^vix'), will play a useful role in the calculation of the coordinate space Feynman diagrams. 
The action (2.1) is invariant under conformal inversions, as [El: 

4>{x) — > 4>'{x) = x' <p(x'), 

ip(x) — > ip'(x) = x' 2 'jrjj;'ip(x'), 
A^x) -> A' pix) = -x ,2 J^(x')A„(x'), (2.6) 
while also the following relations hold, 



d x = — jg- and i = -x' J ^{x')^. (2.7) 

In order to use conformal properties to simplify the two loop Feynman integrals, one 
should expect that the relevant Feynman rules will consist of vertex factors and propagators 



with simple inversion properties. In particular for the scalar and spinor propagators we 
have, 

A(x-y) 



A7T 2 {x-y) 2 4tt 2 (x' -y') 2 ' 
S (X - y) = - mX - y) = = -^M±zi-t. (2.8) 

The vertex rules, read from the action (2.1), are: 




Z2 ' 



where solid lines are fermions, dashed lines are scalars and wavy lines are gauge fields. 

In addition the energy-momentum tensor insertions (2.3) and (2.4) yield the following 
vertices: 

J > z,^v = fc 7( A)° (a| - aff ) 5 \ z ~ *0 <* 4 (* - ^2), (2.12) 

z i x 

7|IV -lh(d a if _a a 1 / a 2 , a 2 \\ 

= ^»M- V — 3 ^ ^a^ 7 ^ 2 u ^ az«a z « 2 ^404 ^ dz»dz v x >> 
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where the double solid lines represent gravitons. 

Let us analyze the conformal properties of the graviton vertices. In order to do that, we 
attach the vertices (2.12) and (2.13) to scalar and spinor legs, and use (2.5), (2.7) and (2.8) 
to obtain, 

S(v-x) 70 A )Q (— - — ) S(x-u) = 

= -x' 8 J-^{x') J v)p {x') { v' 2 j' S(v' - x') 7A (A - A) S (x' - v!) u' 2 i' }, (2.14) 
for the fermionic vertex, where the derivatives only act inside the curly brackets. Likewise, 
A( ) ( ^ ^ ^ ^ d d 1 , <9 2 ^ (9^ ))A( ) J ( J (') 
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- ^ - y„ - 1 ( _ + w - u> *}, (2,5) 

for the bosonic vertex, where once again the derivatives only act inside the curly brackets. 

As an illustration of these coordinate space propagators and vertex rules, we shall now 
look at the one loop triangle diagram and perform the conformal inversion on the amplitude's 
tensor structure. The relevant one loop triangle diagram is depicted in Figure 1(a) and its 
amplitude, B a ^ upa (z,y,x), can be computed using the previous rules to be: 

1 d d d d 

= -igk 2 Tv la ^S{z - y) 7(/ A )/3 (— - — ) S{y-x) ~f ip 5 a)w (— - — ) S(x-z). (2.16) 

Due to translation symmetry we are free to set z = 0, while we refer the remaining external 
points x and y to their inverted images (2.5). Although this transformation may seem ad 
hoc at this stage, it will later simplify the calculation of the two loop diagrams [f|. The 
result we obtain is, 

B a ,iiv, P a(0, V, x) = V 8 J^{y')Jv)u{y')Jp{p{x')J a)5 (x') Tr7 a 7 5 7 P jTfaj S(y'-x') j p . 

(2.17) 

Taking the fermionic trace one finally gets, 



B a)tlu ^a{0,y,x) — ^ g J/ X J[i(p{y ) J u )u{y ) Jp{p{x ) Ja)a{x ) Sqskup Qyi_Q x ,_ ^/ _ ^A4 ' 



(2.18) 



For separated points (2.16) is fully Bose symmetric and conserved on all indices. The 
expected anomaly is a local violation of the conservation Ward identities which arises because 
the differentiation of singular functions is involved || . There are several ways to obtain the 
anomaly in this coordinate space approach [|1(],|14],[| . One way [|iy,fj to do this is to recognize 
that the amplitude (2.16) is too singular at short distances to have a well defined Fourier 
transform. One then regulates, which entails the introduction of several independent mass 
scales. The regulated amplitude is well defined, and one can check the Ward identities. Also, 
an important aspect of this coordinate space approach to the axial anomaly is that the well 
defined amplitude (2.16), for separated points, determines the fact that there is an anomaly 
of specific strength [|l(| . 



Some more comments about the role of the conformal symmetry in the calculation of 
possible radiative corrections to the anomaly are now in order. At first sight this could 
look as a questionable role, after all the introduction of a scale to handle the divergences of 
perturbation theory will spoil any expected conformal properties. This is true in general, 
but our two loop triangle diagrams for this massless Abelian Higgs model are exceptional. 
Any primitively divergent amplitude is exceptional when studied in coordinate space for 
separated points, since the internal integrals converge without regularization ||. 

As we shall see in the next section, there will be 3 non-planar and 3 planar diagrams, 
which are primitives. Of course there will be many other diagrams which contain sub- 
divergent vertex and self-energy corrections, and these require a regularization scale. How- 
ever, we are dealing with pure axial coupling for the fermion. This means that if we introduce 
diagrams containing internal photons, then there is a unique choice of gauge- fixing parameter 
T which makes the one loop self-energy finite 0. Moreover, since the vertex and self-energy 
corrections are related by a Ward identity, each vertex correction is also finite in this same 
gauge. In conclusion, choosing this finite gauge makes it possible to obtain a finite two loop 
result in our calculation. 



S 



III. THE THREE POINT FUNCTION FOR THE TWO LOOP GRAVITATIONAL 

AXIAL ANOMALY 



Let us now proceed to the next loop order in the non-gauge sector, as we are interested 
in computing possible corrections to the gravitational axial anomaly at order 0(gf 2 k 2 ). At 
this order we have a total of 36 diagrams that can possibly contribute. Of these diagrams, 
3 are non-planar, but they actually only correspond to 2 independent calculations due to 
reflection symmetry. These are depicted in Figures 1(b) and 1(c). Then, there are 3 scalar 
self-energy diagrams, and other 3 photon self-energy diagrams (as we shall see, some dia- 
grams involving photons are required in order to choose the finite gauge and compensate 
some divergences of the non-gauge amplitudes). These 6 self-energy diagrams amount to 
2 independent calculations alone, the ones depicted in Figures 1(d) and 1(e). Then, we 
have 3 axial current insertion vertex corrections, in Figures 1(f), 1(g) and 1(h). At the 
energy-momentum tensor insertion, we also have vertex corrections. These are 6 diagrams, 
amounting to the 3 independent calculations in Figures l(i), l(j) and l(k). There are also 6 
diagrams that identically vanish due to fermionic traces, the ones in Figures 1(1), l(m) and 
l(n). Associated to the mentioned self-energies there are 3 diagrams corresponding to local 
self-energy renormalizations. They amount to 1 independent calculation, Figure l(o). Also, 
associated to the mentioned vertex corrections at the energy-momentum insertion there are 
2 diagrams corresponding to local vertex renormalizations. They amount to 1 independent 
calculation, Figure l(p). So, overall, of the 29 initial two loop diagrams in Figure 1, we are 
left with 12 independent calculations. In Figure 2 we have 7 more diagrams, corresponding 
to 4 independent calculations. These diagrams are associated to the finite gauge photons 
and shall be discussed later. We are thus left with an overall number of 16 independent 
calculations, out of the initial 36 diagrams. Let us see how to perform such calculations, 
one at a time. 
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A. Diagrams in Figures 1(b) and 1(c) 



We begin with the non-planar diagram depicted in Figure 1(b), which we shall denote by 
^aJiu, P a( z J Vi x )- This amplitude is conformal covariant since no issues of sub divergences and 
gauge choice arise. The idea M is to use the inversion, u a = u' a /u' 2 and v a = v' a /v' 2 , as a 
change of variables in the internal integrals. In order to use the simple conformal properties 
of the propagators (2.8) we must also refer the external points to their inverted images (2.5), 
as was done in (2.14), (2.15), and in (2.17), (2.18). If in succession we use the translation 
symmetry to place one point at the origin, say z — 0, then the propagators attached to that 
point drop out of the integral, because the inverted point is now at oo, and the integrals 
simplify. 

After summing over both directions of Higgs field propagation, and setting z = 0, the 
amplitude for Nj^ po .(0, y, x) is written as, 

■Tr l5 S(v - x)7(A)/3 (-^r ~ S ( x - U ) S ( U ~ 2/h(/A> - ^0 S(y - v). (3.1) 

The change of variables previously outlined can be performed with the help of (2.7), (2.8), 
and the Higgs current transformation, 

- V U [U a — V a ). (3.2) 



U A V 2 V 4 U 2 



The spinor propagator side factors ^ , j) , etc., collapse within the trace, and the Jacobian 
(u'v')~ 8 cancels with factors in the numerator. Performing the algebra we obtain, 



Ni%J0, y, x) = - l l^-y'\» J a „( V ') J^(y') J f{p (x') J„ )s (x') f dV dV (,/ - 



> 7, SW - x> hf (±- - ±-) Six' - u>)S ( u' - y'H ( £ - j£> S(y' - ,')}, (3.3) 

where the derivatives are acting only inside the curly brackets. 

We see that we obtain the expected transformation factors for the energy-momentum 
tensors at x and y times an integral in which u' and v' each appear in only two denominators. 
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These convolution integrals can be done in several different ways. The final relevant formulas 
are listed in the Appendix. We begin by using the trace properties to move the S(y' — v') 
propagator in (3.3) close to the S(v' — x') propagator. The differential operators are kept 
fixed, with the understanding that now the y' derivative that seems to be acting on nothing 
is actually acting on the propagator S(y' — v') which is now sitting on the left. As usual, 
all derivatives act only inside curly brackets. We can perform the integrations without the 
need to make the differentiations first as the integration variables are well separated from 
the differentiation ones. Expand the product with (v' — u') a , and we are led to the following 
result: 

dV^- M %Tr 75 S(*/-*U 
-Tr 75 / dyv' a S{v'-y')S{v'-x') lp {^-^r) J d*u> 'S(u> -x')S(u> -,fl 7A (JL- ±-) 



+Tr 75 / d A v'S{v'-y')S{v'-x') lp {—-—) J d 4 u'u' a S(u' -x')S(u' -y') lfl ( 



Qy'v dy' P 
(3.4) 



where the integrals can be directly read off from the Appendix. When these results are used 
and substituted within the trace, one finds the final amplitude, 

N al„,pa(0, y, X) = g y' S x' S Jfifaiy') Ju)u{y') Jp(A X ') J a)a(x')- 



32vr 

d d (x' — y') K d' d 



r)}« ( 3 - 5 ) 



(x' — y') 2 dx'a dx'v (x' — y') 2 dy' p dy' 
where we should have the "trace" attitude for taking derivatives: the derivatives only act 

inside the curly brackets, and the y' derivative that seems to be acting on nothing is actually 

acting on the first {x' — y') term. 

We observe that unlike the case in 0, this non-planar amplitude is not a numerical 

multiple of the amplitude for the one loop triangle diagram (2.18). This is because the 

tensor (derivative) structure in (3.5) is different from the one in (2.18). To see that one just 

has to explicitly compute both structures, and compare them. For the triangle one has: 

d 2 A K 
dxLdy' p A 4 
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AAAfii _j_ a. d 2 , i x i a 2 /A„n A rs 6) 

~ ftr^A^d^A 2 / A^x^va 2 ^ AW^^A 2 ' + dx'^ A 2 hy' p ^ A 2 ^ 1 ' J 
while for the non-planar structure one obtains, 

J_ ( J[ t A K t g 

A 2 W ff A 2 l <9y' P 9y'/ 

= Ar^Ar ±_) M.IJLAu^f (37) 

&r'- K -A 2 )dy' u VA 2 ^ A 2 dx'^dyi ^A 2 ' A 2 dx'- a dy'^ A 2 ) dx'^A 2 ) dy'^A 2 )' 1 ' ; 
where we defined A = (x' — y'). The reason such difference can happen is that while in || 
there is a unique conformal tensor structure for the correlator of three axial vector currents, 
in here we have two conformal tensor structures due to the higher dimensionality of the 
correlator of the one axial vector current and the two energy-momentum tensors. Also, 
observe that both these structures (3.6) and (3.7) are to be understood as always attached 
to the appropriate factors of J^ v {y'), J pa {x') and the appropriate powers of y', x'. Moreover 
the diagrams that give rise to them obey conservation equations for the energy-momentum 
tensor insertions. (3.6) is associated to the one loop diagram in (2.18). It can easily be 
proved that the conservation equation is obeyed, a standard result from (and also form 
[pTOfl once we are aware of the relation (A. 6) from the Appendix). (3.7) is associated to the 
two loop diagram in (3.5), and one can also explicitly check the conservation law for this 
case. This existence of two conformal structures is an extra feature in the discussion of these 
two loop diagrams, relative to the work in ||. 

There are 2 more non-planar diagrams, where the scalar vertex is placed at x and at y. 
We need to compute them, as they are independent of the previous result (we have a scalar- 
scalar-tensor vertex instead of a scalar-scalar- vector vertex, among other different vertices), 
but they amount to 1 independent calculation. 

So, we proceed with the non-planar diagram in Figure 1(c), denoted in the following by 
Nf® (z,y,x). The method of calculation is very similar to the one for the previous dia- 
gram, and so we shall perform it in here with somewhat less details. After summing over both 
directions of Higgs field propagation, and setting z = 0, the amplitude for (0, y, x) is 
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written as, 

Jv&^o, v, *) = / « 4 « ^ ^ $^ s(« - *h ( A„ (J - J) s<* - «)■ 

• A( " " v) " 5 ~ 5 ( ^ + ajk» A(!/ - (3 ' 8) 

Performing the conformal inversion is now no harder than it was for the previous diagram. 
The procedure is essentially the same, and if we carry out the algebra we obtain, 

Ni%, p A0, y, x) = y'V 8 J, {lx {y>) J v)v {y>) J p{p (x') J a) ,(x')- 

• J dV dV { Tr 757. S(v' - x') lp ( JjL - JjL) _ „/). 

' A(u " s » " ^ - 5 ( ^ + ^» AW - V) }. (3.9) 

Once again the expected structure emerges, and all we have to do is to perform the integra- 
tions. Using the relevant formulas from the Appendix, we find the final result as, 

N alv,pa(°iy, X ) = ^ 8 4^ 8 y' 8x ' 8 JfiUy') M^V') J P(p( X ') Ja)a{x'){e atpK - 

(x'-y')i , g g A x '-y') K ( t t i t ~S i _g g ^ 

' (z' - y') 2 W ff a*V (z' - y') 2 [ dy\ p dy' p) 2 ^ dy>« dy' n 2 [ dy' p dy' p dy' p dy>J'^ 

(3.10) 

However, one should note the following. In (3.3) both differential operators were first order 
in the derivatives, but in (3.9) the differential operator associated to the vertex (2.13) is 
actually second order. As we expect to have at the end a result similar to (3.5) or (2.18), 
we have to perform one of the derivatives in order for both differential operators to become 
first order. Manipulating this result through a somewhat lengthy calculation, one finds: 

Ni%, pa (0, y, x) = - % j£^ y' 8 x' 8 J-^{y>) J v)v {y>) J p{p {x>) J c) - a {x')- 

J £ 1 ( 1 H A x '-y')^ { !_ 1L)\ ( 311 ) 

where the notation is like in (3.5). One realizes that the structure here obtained is the 
same as in (3.5). So, the 3 non-planar diagrams have the same tensor structure, which is 
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different from the one associated to the one loop triangle. From this result we immediately 
read the last non-planar diagram, the one with the vertex involving the scalar fields and 
energy-momentum tensor located at x. All we have to do is to exchange x with y and \iv 
with pa in (3.11). This actually does not change the amplitude (3.11), so that this third 
diagram contributes with the same amount as its reflection symmetric diagram. 

Finally, we can add these 3 diagrams, and obtain the non-planar contribution to the two 
loop correlator. The overall contribution is simply: 

3 3igf 2 k 2 8 8 

N a ^ U)Pa (0,y,x) = J2 N S^, P a(°,y, x ) = — 64^8~ y ' 8x ' J A(m<V) J »)* (v') j pU x ') J °)v( x 'y 
J £ \ (A *d fro f312l 

B. Diagrams in Figures 1(d), 1(e) and l(o) 

We now proceed to the self-energy diagrams. These will be the same as in the three 
gauge current case ||. We shall see the finite gauge mechanism for the one loop self- 
energies and vertex corrections coming about, as it handles certain divergences by choosing 
a gauge where they are zero |4]||. For this cancelation of divergences we have introduced 
the Abelian field which can be decoupled at the end by setting its coupling to zero. Let 
us see how all that works, by starting with the Higgs self-energy diagram in Figure 1(d), 
and the photon self-energy diagram in Figure 1(e). These are 3 diagrams in Figure 1(d) 
(as we can place the self-energy loop at any of the 3 sides of the triangle), which amount 
to 1 independent calculation, and other 3 diagrams in Figure 1(e) that again amount to 1 
independent calculation. If we remove the self-energy leg from the triangle diagram, and 
add the Higgs and photon contributions we obtain 

E(« - u) = JL If + I/(i _ r)] + a$5\v - it) , (3.13) 

where T is the gauge fixing parameter coming from the photon propagator In this 

result, the first term is the part of the amplitude which is determined by the Feynman rules 
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read from the diagrams. It has a linearly divergent Fourier transform, but the crucial point 
is that this amplitude can be made finite by choosing the gauge T = 1 + 2f 2 /g 2 . It then 
vanishes for separated points. However, there is a possible local term, the second term in 
(3.13), which is left ambiguous by the Feynman rules, and is represented in Figure l(o). The 
constant a will be determined by the Ward identity ||. 

In order to proceed with the calculation of this constant using the Ward identity, we 
first need to look at the following vertex correction diagrams at the axial current insertion: 
Figure 1(f), T^l v>pa (z,y,x), Figure 1(g), Tg^z, y, x), and Figure 1(h), T^l v>pa (z,y,x). 
Again, we need a diagram involving photons in order to choose the previously introduced 
finite gauge. Also, these 3 diagrams clearly correspond to 3 distinct calculations. 

The amplitudes of the 3 vertex correction subgraphs in these diagrams are the same as 
in ||. Therefore we already know that each contribution has a logarithmic divergent Fourier 
transform, and that the sum of the divergent contributions from these 3 vertex subgraphs 
is proportional to —If 2 — g 2 (l — T), therefore vanishing in the same gauge that makes the 
self-energy finite. Henceforth we shall use this gauge. 

Let us then proceed with the Ward identity calculation, by summarizing the result from 
0. From the amplitudes for the vertex subgraphs in the diagrams vpa (z, y, x), i = 1, 2, 3, 
we obtain the Ward identity for the theory J5|, 

d 1 

— T a (z, u, v) = -i-srrs (5 4 (z - u) - 5 4 (z - v)) E(« - v), (3.14) 

where T a = Yh=i T^, and T % a is the vertex subgraph in the diagram T$ . The constant a 
in the self-energy (3.13) can be calculated as in || - where basically one works out the LHS 
in (3.14) (in the finite gauge) in order to find the correct value for (3.13) in the RHS -, and 
the final answer is given by 

nz) = ^ 2 U 2 -\g 2 )^\z)- (3.15) 

Strictly speaking, one should now proceed to verify that the exact same result is obtained 
from the Ward identity associated with the vertex correction diagrams at the energy-momen- 
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turn tensor insertions, Figures l(i), l(j) and l(k). This is in fact true, but for pedagogical 
reasons we shall postpone such a proof for a couple of pages. 

It is this result for E(i> — u) which is to be used to evaluate the local self-energy renor- 
malization, Figure l(o), therefore yielding the correct value for a in (3.13). These again are 
3 diagrams that amount to 1 independent calculation as in Figures 1(d) and 1(e). As (3.15) 
is purely local, the integral in u and v required for the previous diagram is trivial, simply 
yielding a multiple of the one loop triangle amplitude. The final result is that the sum of 
the self-energy insertion diagrams, Figures 1(d), 1(e) and l(o), is a multiple of the one loop 
amplitude, 



3 1 

K,^, P A z ,y, x ) = if - ^9 2 )B a ^ U:P(7 (z,y,x), (3.16) 



exactly like in ]|] as the internal fields are the same. Now recall that there is a factor of 3 
from the triangular symmetry. There is also a factor of 2 for opposite directions of fermion 
charge flow (such term was absent in the non-planar diagrams). Finally, we are interested in 
the 0(gf 2 k 2 ) corrections, so that the term in g 2 in (3.16) should be discarded. The overall 
result for the self-energy contribution to the two loop correlator is finally, 



2 



9/ 

Sq!,^,p<t(0, y, x) = -^—^B a ^ UtPa ( y 0,y,x), (3-17) 
where we have set z = (for coherence with the other diagram calculations). 

C. Diagrams in Figures 1(f), 1(g) and 1(h) 

We can now proceed the calculation of the vertex correction diagrams at the axial current 
insertion, T£> (z,y,x), i = 1,2,3. As for the self-energy, the calculation of these 3 
diagrams follows from ||. We shall regard each virtual photon diagram as the sum of two 
graphs, one with the photon propagator in the Landau gauge T = 1, and the second with 
inversion covariant pure gauge propagator, 

1 f 2 1 

A^iu-v) = — - ^J^iu-v). (3.18) 

47r z g £ [u — v) z 
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The Landau gauge diagrams give order 0(g 3 k 2 ) contributions to the two loop correlator, 
while the remainder gives an order 0(gf 2 k 2 ) contribution which is what we are interested 
in. Therefore - and similarly to what was done from (3.16) to (3.17) - we shall discard the 
Landau gauge diagrams from our final result, and only use (3.18) for the virtual photon 
propagator in the finite gauge. 

With this in mind we turn to the calculation of the diagrams (z, y, x). The method 

of calculation is similar to the one used for the non-planar diagrams, and so we shall follow it 
here without giving details. After summing over both directions of Higgs field propagation, 
and setting z — 0, the amplitude for T^J^, (0, y, x) is written as, 

Tr * 7 " 75 ^ ^ ~ ^<A)/> { dy~, " S{V ' X)l{A> " dx~J ~ M) ' (3 ' 19) 
and performing the conformal inversion we are led to the result, 

n%, p Ao, v, x) = y'V JfUv') MM) JxAz') -M*') / ' <<Va(i/ - »')• 
sw - y'H (|: - |-) W - *h, (|t - |r) - «'>• (3.20) 

As usual the expected tensorial structure emerges. We are left with the integrations to be 
performed. However, as we have seen, this result is divergent; only when we sum the 3 
diagrams Tf£ (z , y , x) , % — 1,2,3 the result will be finite, in the finite gauge. So at this 
stage we should include in the calculation (3.20) the equivalent results coming from the 
diagrams ^(0, y, x) and T^fy (0, y, x) - where for this last one we should use only the 
inversion covariant pure gauge propagator (3.18). The result of including the 3 diagrams all 
together is to produce an integral of a traceless tensor, which is convergent, and can be read 
from the formulas in the Appendix. Hence we can write for the net sum of vertex insertions 
at point z, i.e., T^ pa (0,y,x) plus T^ p(T (0, y, x) plus T^ pa (0, y, x), 

T a,nv,pa(Qiyi X ) = 2567T 8 V % ^ Jv)p(V ) Jp(p( X ) J<r)*( x ) £ aMp.p Qyl Q x l ^/ _ yiy ' 

(3.21) 
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which is a multiple of the triangle one loop amplitude (2.18). Recalling that there is a factor 
of 2 for opposite directions of fermion charge flow, we can finally write for the contribution 
of the vertex correction diagrams (at the axial current insertion) to the two loop correlator, 

P 

(0 iDi x ) — 00 o ^a^u,pa(^,y,x). (3.22) 

D. Diagrams in Figures 2(a) and 2(b) 

In order to obtain the previous result we had to use finite gauge virtual photons, as given 
by the propagator (3.18). If one has a diagram with an internal virtual photon, one should 
expect a factor of g from each of the two internal vertices, and so an overall contribution of 
order 0(g 3 k 2 ). However, if one uses an internal finite gauge virtual photon, there is an extra 
factor of f 2 /g 2 from the propagator (3.18), and we therefore obtain an overall contribution 
of order 0(gf 2 k 2 ), which is the order we are interested in. This means that one now has to 
include all diagrams with one internal finite gauge photon (3.18). 

In particular we have to include one more vertex in our rules, that completes (2.9-13). 
This vertex can be read from the action (2.1) when coupled to gravity, and is the following, 

z i 



vA/W > = Z,|av = -\igk (7 (/ A )a - \ 5^ a ) 75 5 4 (z - z x ) 5 4 (z - z 2 ) 5 4 (z - z 3 ), 
Z, / (3.23) 



where the notation is as in (2.9-13). Observe that there is no similar vertex involving two 
scalar legs (as opposed to the two fermion legs we have) as such vertex would give diagrams 
that do not contribute to the abnormal parity part of the correlator we are computing. 

When this vertex is considered, one finds that there are 7 new diagrams that must be 
included in our calculation, the ones presented in Figure 2. There is 1 primitive diagram in 
Figure 2(a). There are 2 other primitive diagrams which only correspond to 1 independent 
calculation, the one depicted in Figure 2(b). Then we have energy-momentum insertion 
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vertex corrections. These are Figure 2(c) and Figure 2(d), 2 independent calculations corre- 
sponding to 4 diagrams due to reflection symmetry. We shall now proceed to evaluate these 
diagrams. 

We start with the primitive diagram in Figure 2(a), P^]} U:P(T (z, y, x). This diagram is 
easily evaluated as it involves no integrations. Recall that we have to use (3.18) alone, 
whenever one encounters a virtual photon. Setting z = and performing the conformal 
inversion, the amplitude Pa^pai® , y x) becomes, 

P Sv,pa(°, V, X ) = - W q{^8 y' 8x ' S J PUV') J V)AV') Jp(p( X> ) J a)a{ X ') £ a K flp ^ _ ^ Jva(x'-y'). 

(3.24) 

Unlike all the preceeding calculations, this amplitude involves no derivatives. This is 
certainly to be expected due to the nature of vertex (3.23). However, one can manipulate 
(3.24) in order to write it as the second derivative of a tensor involving the structures (3.6) 
and (3.7) alone. After some calculations, one can show that (3.24) can be re-written as 
(including the factor of 2 for opposite directions of fermion charge flow): 

p^I, p Ao, y, x) = ^ y'V 8 j-^{y') J v) M) W) J^{x')- 



J £ 1 ( g d (x^y% B d , 

^ _ y/ y K dx ,_ dx ,_) ^, _ y/ y K dyt _ dy/ _)j, {MO) 



and this tensor structure is precisely the same as the one in the non-planar diagrams (3.12). 
This should not come as a surprise as this diagram - like the non-planar ones - is a primitive. 

Let us proceed with the primitive diagram in Figure 2(b), P^^ pa (z,y, x). This diagram 
involves only one integration, therefore being different from the ones we previously calculated 
(which either involved two or none integrations). Upon setting z = the amplitude for 

PL ( %, P Ao,y,x) is, 

PL%, P A^y,x) = -{^J^yJr,(u-yy 

ik ill 1 

•1*7/3 — 7 a75 — (7(A)r - 2 <Wr) S(y-x) 7(A)tt (^- - ^-) S(x-u), (3.26) 
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and performing the conformal inversion one obtains, 

p'^ p M v, x) = y' 8 *' 8 JM MM) W) J^{x'). 

■ J jj-^y Mu' - V') Tr lp la 75 7 , S(y' - x') lp {±- - ±-) S(x> - u>). (3.27) 

We are left with one integration to perform. However, one should note that there is only 
one differential operator in (3.27), and if we are to obtain a final result for this diagram which 
involves the tensor structures (3.6) and (3.7) we shall have to manipulate (3.27) in order to 
re-write it in such a way that it involves two differential operators. This is analogous to the 
situation we faced from (3.24) to (3.25). After integrating and performing some calculations, 
one obtains: 

P a%,pa(°, V, X ) = t J^8 y'^' 8 J pUy') J u)u{y') Jp(p(x') Ja)a( X ') e <*WP 

r 1 d 2 (x> - y_% 1 ( J[ (x' - y% ][_ _ ^d_ , 

^2 dy' v dx' & {x'-y'Y {x 1 - y') 2 V dx' 9 dx' 5 ' (x 1 - y') 2 { dy' P dy'J* K ' ' 

It is interesting to observe that this tensor structure is a linear combination of both (3.6) 
and (3.7). Also, one must now include in this result a factor of 2 for opposite directions 
of fermion charge flow and another factor of 2 associated to the 2 distinctive diagrams 
connected through reflection symmetry. 

Finally, we can add the 3 diagrams in Figures 2(a) and 2(b), in order to obtain the 
primitive planar contribution to the two loop correlator, 

igpfc 2 g g 

Pa,pu,piri^i y-i x ) = 22^8 y X Jp-faiy ) Jf)v{y ) Jp(p( x ) Ja)a{% ) ^anjip" 

rl d 2 (x'-y') K 5 1 g t (x'-y>) t t } 

b dy' D dx'- a {x'-y'Y 4 (x' - y') 2 K dx'„ dx' J {x' - y') 2 K dy' p dy'J> K ' ' 

E. Diagrams in Figures l(i), l(j), l(k), l(p), 2(c) and 2(d) 

Next, we proceed with the evaluation of the contributions coming from the vertex correc- 
tion diagrams at the energy-momentum tensor insertions. These amplitudes are presented in 



20 



Figure l(i): V^ Utf)tr (z, y, x), Figure l(j): V^ pa {z,y,x), Figure l(k): V^ pa (z,y,x), Figure 

!(P) : v l%, P a ( z , V, x ), and also Fi S ure 2 ( c ) : v l%, P a( z , V, x ), and Fi S ure 2 ( d ) : V a%, P a( z , V, x )- 
Using the previous treatment with the axial insertion vertex we should insert the photon 

propagator in our calculation in order to guarantee finiteness of the energy-momentum in- 
sertion vertex. The amplitudes and have a different structure from the other vertex 
diagrams as they are semi-local. Moreover, there is a possible local term, which is left 
ambiguous by the Feynman rules. This is analogous to the situation we faced when dealing 
with the self-energy diagrams. As this local term cannot be evaluated by Feynman rules it 
will be solely determined from the Ward identity. The role it plays is one of regularizing a 
divergence. 

The sum of these amplitudes becomes finite in the finite gauge, which also guaranteed 
the finiteness of the axial insertion vertex. The finite part of the vertex subgraphs is a 
traceless tensor with respect to all three included indices, and can be written as, 

v^(v v u) = 3k(f2 - 1 7 ( (y-u),(y-u)Ay-u) K 1 x 

(3.30) 

In addition, it would be thoughtful to check the Ward identity connected to these vertices. 
For that we express the tensor on the RHS of (3.30) in terms of the regularized traceless 
structure of derivatives, 

( {y-u)fi{y-u) v {y-u) K 1 ( , \ _ _ LJLfJLJL I a n\ 1 

^ (y - uf 6 ^ [V UU) ) - 48 dy K ydy, dy v 4 V U ) (y - uf ' 

(3.31) 

With this expression it is easy to derive that the following Ward identity is satisfied for the 
energy-momentum insertion vertex, 

^ J d*vj2V$(v,v,u) = -kd v Z(y - u), (3.32) 

where VjfJ is the vertex subgraph in the diagram . 

This Ward identity is essential for the calculation of the amplitude of Figure l(i). It 
suggests that it will give a divergent result, and only when we add together the diagrams in 
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Figures l(i), l(j), l(k), l(p), 2(c) and 2(d) we shall obtain a finite answer. After summing up 
the two possible directions of the Higgs field and setting z = due to translation invariance 
of the amplitude, V^ upa (0,y,x) is written as, 

igfk 2 



^2^(0,0,*) = J d*ud%A(u-v) 



Tr S(v - y) 1{ ,5 u)p (A - A) S (y - u)S(u - x) 1{p 5^ (A _ JL) 7a7s (3.33) 

Using the conformal properties of the theory we can perform the usual inversion in the 
spatial variables which results to, 

K2^(0, y, x) = ^A y'V 8 J^(y') -M-Av ) JpUx') Ja),(x') I dV dV A(t/ - u')- 

•Tr 75 S{v' - y'h, (A _ JL) S {y> - u')S(u> - x') lp A 7a . (3.34) 

In order to proceed with the if and v integrations we should also add the diagrams of Figures 
l(j), l(k), l(p), 2(c) and 2(d) to guarantee finiteness of the contribution from the energy- 
momentum insertion vertex. Note that all these diagrams do not contribute with a finite 
part for the order 0(gf 2 k 2 ) we are interested in, but merely make the diagram l(i) finite. 
This will make the integrand have a traceless form, as given in the Appendix. After some 
manipulations we deduce that, 

Va„& Vi X) = y' 8 x > 8 J_ (/i ( y ') J u)p (y') J_ (p ( x ') J(j) -( x >) E(XKpp - AA |^ _ ^ , 

(3.35) 

which is proportional to the triangle structure. Taking into account the 2 fermionic directions 
and doubling our answer for the two distinctive diagrams connected with reflection symmetry 
we obtain finally: 

2f 2 

(0, y, x) = -^-^ B a ^ UtPa (0, y, x), (3.36) 

which is similar to what we got in (3.22). So, we can add all the diagrams that represent 
vertex corrections (both at axial and energy-momentum insertions). The overall result of 
the vertex corrections contribution to the two loop correlator is, 

3f 2 

(0 1 x ) — 00 2 ^a,fiu, P cr{0, y, x). (3.37) 
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F. Diagrams in Figures 1(1), l(m) and l(n) 



Finally, we would like to mention the diagrams that are zero. That these diagrams vanish 
can be easily seen either from the fact that the fermion trace vanishes or from arguments of 
Lorentz symmetry. We mention these diagrams for completeness. They are the following: 
Figure 1(1), which are 3 diagrams that amount to 1 independent calculation, Figure l(m), 
which is only 1 diagram, and Figure l(n), which are 2 diagrams that amount to 1 independent 
calculation. We have now completed the calculations for all the 36 diagrams. 



G. The Three Point Function 

The next and final step is to add all diagrams together, and find out what is the two loop 
contribution to the three point function at order 0(gf 2 k 2 ). Adding the results for all our di- 
agrams we obtain the 0(gf 2 k 2 ) two loop contribution to the correlator (A a (z)T^ u (y)T pa (x)) . 
There are 4 distinct contributions: the one from the non-planar primitive diagrams, N atmpa 
(0,y,x) in (3.12); the one from the self-energy diagrams, £ ajAO/)Pcr (0, y, x) in (3.17); the one 
from the planar primitive diagrams, P a ^ U)Pcr (0,y,x) in (3.29); and the one from the vertex 
correction diagrams, V atPU ^ pa (0,y,x) in (3.37). Adding these 4 structures we finally obtain 
our result: the three point function does not vanish and consists of two independent confor- 
mal tensor structures, 



igpk 2 

64vr 8 



y' S x' S Jp{p{y') Ju)u(y') Jp{p{x') Ja)a(x') 



'ClKpp 



f 7 d 2 (x'-y') K 11 1 g t (x'-y') t 

L dy'^dxL (x 1 - y'Y 2 (x 1 - y') 2 W s dx' s ' {x< - y') 2[ dy' D dy'J* l ' ' 

There is one consistency check that can be performed on this result. Namely, under the 
appropriate changes, one can ask: does it reduce to the result obtained in the axial gauge 
theory case 0? In order to reduce (3.38) to the gauge theory case of we first have to 
discard all diagrams involving the vertex (3.23). Then, in the other diagrams, one has to 
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erase all the "graviton derivatives" . Once this is done we are left with a unique conformal 
tensor, i.e., we are reducing the structure of our theory to the one in ||. Finally, taking into 
consideration that the removal of the "graviton derivatives" also includes a factor of —2 in 
the non-planar contribution (due to the symmetry enhancement of this diagrams), one finds 
that the overall result vanishes just as it did in |J. This shows that our result is consistent 
with the calculations performed for the gauge axial anomaly. 

We can trace back the reason why this radiative correction does not vanish (as it does 
vanish in the gauge theory case ||). This is (3.12) and (3.29), the contribution of the 
primitive diagrams, which is not a multiple of the one loop amplitude. The existence of two 
different conformal tensors in our theory is a result of the dimensionality of the correlator 
(A a (z)T^(y)T p(7 (x)). 
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APPENDIX A: DIFFERENTIAL REGULARIZATION 



In this Appendix we study the differential regularization of the one loop triangle diagram 

associated to the gravitational axial anomaly, Figure 1(a). As we have seen in Section II. 

the amplitude for this diagram is, 

1 d d d d 

-igk 2 Tr 7a75 S(z - y) 1{ll S v)i (—-—)S(y- x) 1{p 5 a)j (— - — ) S(x - z). (A.l) 

We can re-write this diagram as a "generic" fermion triangle diagram, for which the bare 
amplitude takes the form: 

s{z - i - h s{v ~ x)1 - i ~ i? s(x ~ z) - (A - 2) 

As compared to (A.l), this is a slightly more general form of the amplitude we want to 
regulate. For our present purposes we shall only need to concentrate on the regularization 
of the singular functions present in the amplitude, so we may as well just consider (A. 2). 
Performing the derivatives, (A. 2) can be written as: 

d 2 

Tr { -q^t_ (im S(z - y) 7// S(y - x) 7/ S(x - z)) + 

d d 
+ 27/j — S(y-x)^ I — S(x-z) j in S(z - y) + 

d d 
+ 27m -S(z - y) 777 — S(y - x) 77 S(x - z)} = 

1 r d 2 / d 1 d 1 d 1 \) 

~ "(4^2)3 [7/// la III lb II 7cJ ( -fry-faTX dz~ a {^ z - y)2 dy~ b (y - X f di~ c ( X - zf ) J 

2 r g g 1 g g 1 g 1 . 

-(4^2)3 [77/ 7a 7/ 7b 7m 7cj 1 7^ 7^ 7^—^ 7^ 7^ 7^77-^ 7^ ) ~ 

2 r g g 1 g g 1 g 1 1 

- (4^2)3 iY 7a 7// 76 li 7cJ ( t^T ^" ( z _ y )2 3^ ( y _ x )2 7^ ( x _ ^2 /• ^) 
With x' = x — z, y' = y — z, a = c,b = a and c = b, we can re-write the first term in the 
previous result as the singular function, 



(1) / / a d 2 / d 1 d I d 1 \ 



though in the following we shall drop primes and bars, and simply use the notation (x, y) 
and {a,b,c}. By a similar re-naming of variables one can likewise manipulate the second 
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and third terms in (A. 3) to obtain - in both cases - the singular function, 



l 3 abc 



d d 1 d d 1 d 



dxj dx a x 2 dy { dy b y 2 dx c (x - y) 2 ' 



(A.5) 



Observe that (A.5) is independent of (A. 4), so that in (A. 3) we have two independent 
singular functions. For particular choices of the vertex gamma matrices, (A. 2) describes the 
anomalous three point function (explicitly studied in the paper) as well as other physically 
interesting amplitudes. 

We need to "pull out" derivatives and regularize both ,{x, y) and . (x, y). Start- 
ing with (A. 4), one easily observes that this can be written as, 

d 2 



t 



% 3 abc 



x,y) 



tabc(x,y) } 



(A.6) 



dyidxj 

where t ab c(x, y) is the singular function which is present in the bare amplitude for the gauge 
axial anomaly, and has been previously considered in |H| . In particular, it was shown in this 
reference how to regularize this singular function. Two derivatives are required to control 
the linear divergence arising from the singularity at x ~ y ~ 0. Making manifest the x <-> y, 
a <-» b antisymmetry of t a b c (x,y) one can write, 



t a bc(x,y) = F abc (x,y) + S abc (x,y), 



(A.7) 



where the function F abc (x,y) has finite Fourier transform by power counting and trace ar- 
guments, 



d d r 1 d 
F abc (x,y) = — - — 



+ 







i_r. 



d 2 



- -S br n 



dx a dy b L x 2 y 2 dx c (x — y) 2 J dx a L x 2 y 2 ^ dx b dx c 4 



d_ r J_ / d 2 _ 1 p x 1 

dy b L x 2 y 2 ^ dx a dx c 4 ac > [x — y) 2 



1 



(x - y) 5 



x 2 y 2 



9 3 



l fs d . d 
- ( o ab - h d bc - 



r 5 , 
OX) 



(x - yf 



(A.8) 



ox a ox b ox c o v aa; c ax a 
The term S abc (x,y) contains the "traces" subtracted off in (A.8) and thus derivatives of 
5(x — y) times 1/x 4 factors which are regulated as is standard in differential regularization, 
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yielding: 



,v) = -y{ 



>bc 



_d_ 

dx n 



- 5n 



d_ 

oy b 



5(x -y)U 



In Mix 2 



or 



Obc 



OT a oy a ' v dx 6 c% ' v dx c dy c > J x' 
Two different mass scales were used for the two independent trace terms in S a b c (x,y). 
Renormalization or symmetry conditions may be used to determine the ratio Mi/M 2 in 
particular cases of the triangle amplitude. 

Expressions (A. 8) and (A. 9) provide the required regularization of (A. 4) via expressions 
(A. 6- 7). We are thus left with the regularization of (A. 5), which can be performed in a 
similar fashion. Again, one can write, 



d d N1 „, s _ mMfx 2 



)\S(x-y)D 1— }. (A.9) 



(A.10) 



The function {x, y) whose Fourier transform is finite by power counting and trace 
arguments is, 



d 



, . d d d d r 1 <9 

ijahA x i y> ~ "nTT a™ a™ fl„. 



1 



dyi dxj dx a dyb ^ x 2 y 2 <9x c (x — y) 2 J ch/j <9x 



+ 



9 d 



I ^ o 2 



l 



dx a L x 2 y 2 ^ dxbdx c 4 
1 r «9 3 



1 



1 



(x — y) 2 -! dyb*- x 2 y 2 ^ dx a dx c 4 
(9 .9 .9 



(x - y f 



x 2 y 2 


- dx a dxbdx c 




d d , 
% <9x a ' 




1 


r <9 3 


x 2 y 2 


- dxbdxjdx c 
d d , 


1 


r <9 3 



1 / . a .a . a \ n i 1 -i 

a ( OabT. \- ObcT. V o ac — ) □ — \ + 

h \ r/T_ OXh '' r - — " ' 



5x a 

1 / <9 2 



(x - 



x - yf 



1 / . d d d 



dx h 



dx-i 



(x - yf 



•2„2 V 



<9 2 



x 2 y 2 



dyb dxj <9x a L x 2 y 2 ^ dxidx c 4 ' (x — |/) 2 

1 / „ . „ . „ 

<9x 



ox a oxiOx c o v <7X C 



- ( '\/. — h 4 — 1" 4, — ) D 



{4- 



d 3 



_d__d_ 

dyi dy b 1 x 2 y 2 

d d r 1 

dxj dx a x 2 y 2 L dxidxbdx c 6 ^ to <9x 



dx a dxjdx c 6 
<9 3 



5 r 5 

+ ^cT^ h Q 



x - y)' 



dx c 
d 



dx a 

+ ^6c q h A 



c 9xj 
g 

<9xft 



□ 



□ 



}- 

_L | 

(x - y) 2 J 

-7/^2 } ~ 



(x -y) 2 



27 



o o i 1 



dx a dy b L x 2 y 2 



<9 3 



dxidxjdx c 



1 / d d d 



<9av 



□ 



+ Terms invovling 4 th and 5 th order traceless derivatives acting on 



(x - y) 2 
1 



} + 



(A.ll) 



(x - y) 2 

The fourth and fifth derivatives can be obtained in a straightforward fashion, but their ex- 
plicit form is not relevant and they would occupy a couple of pages to write down. Therefore 
we omit these terms. 

Again, the term S' ijabc (x,y) contains "traces" subtracted from (A.ll) and so its structure 
is similar to the one of (A. 9), containing the usual differential regulated derivatives of 5(x— y) 
times 1/a; 4 factors. One obtains, 



_d__d_l 2 f 

dyi dxj 4 l 



r 9 s 9 1 u ,\nM 2 x 2 
5bc dx- a - 5ac dy-J 5{X ' y) 



f he ( 



<9a; a dy a 



) + 5 ac ( 



(9a; 



d d 1 * r ~ d , _ lnAffx 2 



<9y 6 ' ab ^ dx c 
d 



d \ i n/ x _ \nM 2 x 2 



dy c 



\1 r , , ,_, 111 JWoX -I 



% <9a; a 4 L dy b 



+ 



d 



(9a;/, 



) 



+ 



<9a; r 



9 N1 „, , _ lnM 4 V 



dy c 



)]S(x-y)D 



x' 



! 9 9 1 2 r . 9 . mMfx 2 



9a; n 



* c ^ (9a; a <9ya ^ ac ^ '^' r ' ^ a * ^ '^' r - < ^ ,/ - ^ ~ ^ ^ 



9 9 7T 2 



+ 



% dy b 12 

9 9 7T 2 



(9 



5 



dyi 
d 



dy c 



x* 
In Mix 2 



x/ 



} 



9 



9 



d 



U U TV f r / U U \ / U U \ / U U \"| . 



cfe a 12 



~dx a dy a 

{ K - w) +6lc " as) ^ " ^ "^H 



d d n 2 ( r ^ / d d \ . / 9 d \ . / d d \i . In M|x 2 1 

^ ^ 12 i ^ c W - +5iC ^ " - + 

+ Traces subtracted from (A.ll) in the 4 </l and 5 th order derivatives, 



^<9a; c 

-dx c 
<_d_ 

^dXr 



dy c 
d_ 

dy c 



□ 



In M 2 x 2 



and their respective mass scales. 



(A.12) 



One can see that several different mass scales were introduced for the several independent 
trace terms in Slj (x,y). Again, as was mentioned for (A. 9), renormalization or symmetry 
conditions may be used to determine the ratios between the mass scales in particular cases 
of the triangle amplitude. 
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APPENDIX B: CONVOLUTION INTEGRALS 



In this Appendix we list the convolution integrals that are required in order to perform 
the two loop calculation |15| |. They are inclosed to make this work self contained for the 



reader who wishes to reproduce our result. The table of convolution integrals is (defining 
A = x — y, and using the cutoff A): 

(B.l) 



v 2 (v — x) 2 A 2 



( V X )p jA„. „2 X P 



d% = -n 2 =2 (B.2) 



v 2 (y — x) 4 x 



{v-x) p {v-y) A u 2 A P A mq x 



(v - x) 4 (v - y) 4 2A 2 y pa A 2 



VpVa ~ a5puV2) d*v = £ (x p x„ - -5 pa x% (B.4) 



(v - x) 2 2x 2 v p 4 

(Q - x) p (v - x) a - \5 pa {v - x) 2 ){v - y) x 
(v - xY (v - yY 



(5 pX A a + 5 aX A p - 2^|^), (B.5) 



- - — — - r/V = ^— ir p r n - \>S f ,, : r 2 ). (B.6) 



v G (v - xf 2x 4 v p 4 

((v - x) p (v - x) a - \6pqiy - xf){v - y) x 
(v — x) 6 (v — yY 

-^(^ + ^+1^-4^1^), (B.7) 



d'v 



I X v -*nv-$ * v = i& {6 °> A ° - S - A ' + {x + y) " [s " - 2 ^ D ' (B ' 8) 
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